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Abstract 

o 

' The Bn hyperbolic Sutherland spin model is expressed in terms of a suitable set of 

commuting Dunkl operators. This fact is exploited to derive a complete family of 
commuting integrals of motion of the model, thus establishing its integrability. The 
Dunkl operators are shown to possess a common flag of invariant finite-dimensional 
linear spaces of smooth scalar functions. This implies that the Hamiltonian of the 
\^ • model preserves a corresponding flag of smooth spin functions. The discrete spec- 

trum of the restriction of the Hamiltonian to this spin flag is explicitly computed by 
| triangularization. The integrability of the hyperbolic Sutherland spin chain of Bn 

, type associated with the dynamical model is proved using Polychronakos's "freezing 

Q ■ trick" 
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I. INTRODUCTION 



Since the publication of the pioneering papers of Calogero || and Sutherland j32| |33[ , 
the study of solvable and integrable quantum many-body problems has become a fruitful 
field of research with multiple connections in many branches of contemporary mathematics 
and physics. From a mathematical standpoint, one of the key developments in the field was 
■ the discovery by Olshanetsky and Perelomov of an underlying An root system structure for 

both the Calogero and Sutherland models pj| . The integrability of these models follows by 
expressing the Hamiltonian as one of the radial parts of the Laplace-Beltrami operator in a 
symmetric space associated with the given root system. It was also shown in this paper that 
the original inverse square (Calogero) and trigonometric/hyperbolic (Sutherland) potentials 
arise as appropriate limits of the most general potential in this class, given by the Weierstrass 
p-function, and that integrable models associated to other root systems also exist. The 
rational and trigonometric Calogero-Sutherland (CS) models are also exactly solvable, in 
the sense that their eigenfunctions and eigenvalues can be computed algebraically. In fact, 
the study of the eigenfunctions of these models has led to significant advances in the theory 
of multivariate orthogonal polynomials [Q, [HJ ^3|. Apart from their mathematical interest, 
CS models have found numerous applications in diverse areas of physics such as soliton 
theory J22|, ^9) , fractional statistics and anyons |4|, [| , random matrix theory |35| , and Yang- 
Mills theories || to name only a few. 

During the last decade, CS models with internal degrees of freedom have been actively 
explored by a variety of methods, including the exchange operator formalism [ p4[ , the Dunkl 
operators approach |2|, |sj, reduction by discrete symmetries [^0|, and construction of 
Lax pairs |ll], |2^, pi] . Historically, the first CS models with spin discussed in the literature 
were related to the original models of An type introduced by Calogero and Sutherland |2| 
0, |(], H HI- The integrability of these CS spin models was established in some 
cases by relating the Hamiltonian to a quadratic combination of Dunkl operators of An 



type tj\, 10, 26 1 . The Bn counterpart of the An CS spin models mentioned above were 
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first considered by Yamomoto 38 . In this paper, the spectrum of the rational spin 
model was explicitly determined, and its integrability was shown by means of the Lax 
pair approach. In Ref . J39| , Yamamoto and Tsuchiya presented an alternative proof of the 
integrability of this model using Dunkl operators of -B/v type. The same operators were later 
employed by Dunkl to construct a complete basis of eigenfunctions In contrast, the 

trigonometric/hyperbolic Bjy spin model has received remarkably little attention. In our 



recent paper [|L3| we proved that this model is exactly solvable in the sense of Turbiner 1 36 



|37f| , meaning that its Hamiltonian leaves invariant a known infinite increasing sequence 
(or flag) of finite-dimensional linear spaces of smooth spin functions. In fact, in jl3| we 
developed a systematic method for constructing exactly (or in some cases partially) solvable 
S/v-type CS models with spin by combining several families of Dunkl operators. The key 
elements of this method — first introduced in the A^ case in — are: i) the definition of 
a new family of Dunkl operators, and ii) the construction of a very wide class of quadratic 
combinations of these operators and those in the other two families considered by Dunkl 

The interest on CS spin models has been further enhanced by their close connections 
with integrable spin chains of Haldane-Shastry type |l) . Spin chains describe a fixed 
arrangement of particles that interact through their spins. A well-known example is the 
Heisenberg spin chain, whose spins are equally spaced and interact only with their nearest 
neighbors. The Haldane-Shastry model was actually the first one-dimensional spin chain 
with long range interactions whose spectrum could be computed exactly. In this model, 
the spin sites are equally spaced in a circle and interact with each other with strength 
decreasing as the inverse square of the chord distance between the sites. The integrability of 
the Haldane-Shastry spin chain was proved by Fowler and Minahan in [fblf . Polychronakos 
later realized that the commuting conserved quantities of the Haldane-Shastry spin chain 
can be elegantly deduced from those of the (dynamical) Sutherland spin model of An type 
by applying what he called the "freezing trick" |^7j (see also Q). This corresponds to 
taking the strong coupling limit in the Sutherland spin model and restricting to states with 
no momentum excitations, so that the internal degrees of freedom remain the only relevant 
variables in the problem and the particles are "frozen" at their classical equilibrium positions. 
This observation is, in p rinciple, valid for any integrable spin Calogero-Sutherland model. 
For instance, in Ref. [g7| the freezing trick is applied to the spin Calogero model with rational 
interaction to construct a new integrable spin chain of rational type in which the sites are 
no longer equally spaced. The spectrum of this chain was later calculated by Frahm jl5) 
and Polychronakos p8[ . Bernard, Pasquier and Serban j3j studied the spin chain associated 
with the trigonometric Sutherland model, establishing its integrability for certain values of 
the parameters in the Hamiltonian. 

The aim of this paper is twofold. In the first place, we prove the integrability of the 
hyperbolic Sutherland spin model of Bn type, from which we are also able to deduce the 
integrability of the spin chain associated with this model. Secondly, we give an explicit 
formula for the eigenvalues of the dynamical model whose corresponding (square- integrable) 
eigenfunctions lie in the invariant flag mentioned above. The paper is organized as follows. 
In Section [o| we introduce a family of commuting Dunkl operators of Bn type. We show 
that the sums of even powers of these Dunkl operators generate a complete set of com- 
muting integrals of motion of the Hamiltonian of the model. The commutation relations 
satisfied by the Dunkl operators and the usual permutation and sign reversing operators, 
which possess a richer str ucture than in the rational case, play a key role in the proof of 
this result. In Section HI we analyze the spectrum of the Hamiltonian for any value of the 
spin. Our analysis is based on the fact that the Dunkl operators leave invariant a flag of 
finite-dimensional linear spaces of smooth scalar functions. This flag yields a corresponding 
invariant flag of smooth spin functions for the hyperbolic Bn Sutherland Hamiltonian. We 
construct a partially ordered basis of this "spin" flag in which the Hamiltonian is repre- 
sented by a triangular matrix. In this way we can explicitly compute the eigenvalues of the 
restriction of the Hamiltonian to the (finite-dimensional) intersection of the spin flag with 
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the Hilbert space of the system. We shall use the term algebraic in what follows to refer to 
the these eigenvalues and its corresponding eigenfunctions. It remains an open problem to 
determine whether the algebraic sector of the spectrum actually coincides with the discrete 
spectrum. We also study in detail the (algebraic) ground state, determining its degeneracy 
for all values of the spin. In Section IV we define the spin chain associated with the hyper- 
bolic Sutherland spin model of Bn type, and apply the freezing trick to derive a complete 
family of commuting integrals of motion of this chain. 



II. INTEGRABILITY OF THE SUTHERLAND SPIN MODEL OF B N TYPE 

The Hamiltonian of the hyperbolic Bn Sutherland spin model is defined by 
H* = - d xi + 2 « J2 [ sinh_2 x ij (« + Sij) + sinh" 2 x+ (a + Sij) 

i i<j 

+ b ^2 sinh -2 ^ (b + Si)-b'^2 cosh -2 ^ (b 1 + Si) , 



(1) 



where xh = Xi±Xj and a, b, b' are real parameters. Here and in what follows, any summation 
or product index without an explicit range will be understood to run from 1 to N, unless 
otherwise constrained. The operators Sij and Si in Eq. act on the finite-dimensional 
Hilbert space 

S= (\s u ...,s N ) | Sl = M. M ■ 1 \l: M e ^n), (2) 

associated to the particles' internal degrees of freedom, as follows: 



Sij \s\, . . . , Si, . . . , Sj , . . . , Sn) — |si,...,<Sj',...,Sz,...,Stv), 

SAsi, . . . , Si, . . . , s N ) = |si, . . . , -Si, . . . , sn) ■ 



(3) 



We have also used the customary notation Sij — SiSjSij. 

The operators Sij and Si are represented in S by (2M + l) Ar -dimcnsional Hcrmitian 
matrices, and obey the following algebraic relations: 



S^ — 1, SijSjk — SikSij — SjkSik, SijSki — SkiSij, 
Si = 1 , Si Sj = Sj Si , Sij Sk — Sk Sij , Sij Sj = Si Sij , 



(4) 



where the indices i, j, k, I take distinct values in the range 1, . . . ,N. The algebra 6 generated 
by the operators , Si is thus isomorphic to the group algebra of the Weyl group Wjv of 
type Bn, also known as the hyperoctahedral group. 

We shall also make use of the permutation operators Kij = Kji and the sign reversing 
operators Ki (1 < i ^ j < N), whose action on a function /(x), with x = (x\, . . . , xn) G 
is defined as follows: 

\Kijf) [x\, . . . , Xi, . . . , Xj , . . . , xn) = f (x\, . . . , Xj , . . . , Xi, . . . , xn) , ,g. 
(Kif)(xi, . . . ,Xi, . . . ,x N ) = f(xi, ■ ■ ■ , -Xi, . . . ,x N ) ■ 

The operators and Ki obey algebraic identities analogous to (||). We shall denote by 
& ~ 6 the algebra generated by the coordinate permutation and sign changing operators 
and Ki. Note also that the operators n,j = KijSij and = KiSi generate an algebra 
2B isomorphic to R and ©. From now on we shall identify the abstract group Wn with 
its realizations generated by the operators Kij,Ki on C co (M. N ), Sij, Si on S, or n^Hi on 
C°°(M. N ) (g) S, depending on the context. 
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The Hamiltonian ([!]) describes a system of TV identical particles, whose physical states 
are therefore either totally symmetric or totally antisymmetric under particle exchange. 
Moreover, since H* clearly commutes with the family of commuting operators Hi (i = 
1,...,N), we can choose a basis of common eigenfunctions of H* and all the operators 
11^ . Given an element ipk of this basis, it follows from the commutation relations of the sign 
reversing operators H with the permutation operators Ily that Hiipk = efeV'fci independently 
of i. In principle, the parity could depend on k. However, we shall see in the following 
section that all the algebraic eigenfunctions have the same parity, and that this parity is 
determined by the sign of the parameter 6. From now on we shall assume, for definiteness, 
that we are dealing with a system of fermions whose algebraic states are also antisymmetric 
under sign reversal of each particle's spatial and internal coordinates. This covers what is 
perhaps the physically most interesting case, namely that of a system of spin 1/2 particles, 
for which the internal degrees of freedom are naturally interpreted as the particles' spin. 
The results of this paper can be easily modified to treat any other choice of the particles' 
statistics and parity. 

In the rest of this section we shall prove the integrability of the model (|l|) by expressing 
the Hamiltonian in terms of the following family of B^-type Dunkl operators: 



Ji = dxi - [O- + cotha^iTy + (1 + coth xf- ) K v 



- [6(1 + cotha^) + 6' (1 + tanhx;)]^ + 2a V"ify , (6) 

j<i 

where = KiKjKij. The operators Ji are related to the operators introduced by Ya- 
mamoto J38| in connection with the trigonometric BCn spin Sutherland model. 
A key property of the Dunkl operators (Q) is their commutativity: 

[Ji,Jj]=0, i,j = l,...,N. (7) 

We shall also make use of the following commutation relations between the operators , Ki 
and the operators Ji. 

\ts 7i \2a{K jk - K ik )Kij , i<k<j 

[Kij,Jk\ = < il , (») 

, otherwise 



KijJi - J,K l3 = -2a(l + K H K ii) - (9) 

i<Kj 

KijJj - JiKtj = 2a (l + K v K ji) > ( 10 ) 

i<l<j 

[K h Jj] = 2aK tJ (K 3 - Ki) , [Kj, Ji] = , (11) 

{K h Ji} = -2(6 + 6') ~2aJ2 Ku{K t + K t ) , (12) 

l>i 

where i < j and k are three distinct indices ranging from 1 to JV. Note that the Dunkl 
operators of rational type considered in Ref. |3^] satisfy the latter equations with the right- 
hand side replaced by zero. The nonvanishing of the right-hand side of Eqs. (||)-([l2j) gives 
rise to some nontrivial technical points in the proof of the integrability of the Hamiltonian 
(|l|), as we shall see below. 

Another important property of the Dunkl operators Ji is that they preserve the linear 
space 

Tim = (/x(x) exp [2 njXj) n t = -m, —m + 1, . . . , m, i = l,...,NJ, (13) 
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where 



/i(x) = | sinhx^ sinhxj| a • J^J | sinha^l [coshx^] , (14) 

i<j i 

for any nonnegative integer m. This fact will prove crucial for the calculation of the spectrum 
of the hyperbolic Bm Sutherland spin model. 

We define a mapping *: Q®&^ where 2) denotes the algebra of linear differential 

operators on C ca (M. N ), as follows. If D G D, we set 

(D K ai --- K ar )* = (-l) r D S ar ■ ■ ■ S ai , (15) 

where ctk stands for ij or i. The mapping * is then linearly extended to J) For instance, 
the Hamiltonian of the hyperbolic Bn Sutherland spin model ([j]) is obtained by applying 
the "star" mapping to the operator 



H = - 9l + 2a [ sinh ~ 2 x ij ( a ~ K v) + sinh_2 x tj (« ~ K l3 ) 

i i<j 

+ b^ sinbT 2 ^ (b - Ki) cosh" 2 ^ {b' - K t ) 



(16) 



Let Ao be the antisymmctrisation operator, defined by the relations Aq = Ao and Hy Ao = 
— A , j > i = 1, . . . ,N. More explicitly, 

1 m 

i=i 

where Pi denotes an element of the realization of the symmetric group generated by the 
operators Ely, and e; is the signature of P/. For instance, if N = 2,3 the antisymmetriser 
Ao is given by 

N = 2 : A = i (1 - n 12 ) , 

n = 3 : A = - (i - n 12 - n 13 - n 23 + n 12 n 13 + n 12 n 23 ) . 

The total antisymmetriser A with respect to the action of the operators EEy and IT is 
determined by the relations A 2 = A and 

Ely A = -A, ITA = -A, j > i = 1, . . . , N. (17) 

It may be easily shown that 

A =^(rK i -^)) A o- 

Since Kfj = Kf = 1, the relations ([l7|) are equivalent to 

KijA = -SijA, KiA = -SiA, j > i = 1, . . . , N. (18) 

^From these relations and the definition of the star mapping it follows immediately that 

AA = A*A (19) 

for every operator A G The proof of the integrability of the hyperbolic Bn Sutherland 

spin model is based on the following lemmas. 
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Lemma 1. If B e £> ® S satisfies BA = 0, then B = 0. 
Proof. The operator I? S 5} ® 6 is of the form 

S = ^/ l (x)B i a l , (20) 
iei 

where Bi 6 6, /j 6 C°°(K ), i = (ii, . . . ,ijv) is a multiindex belonging to a finite subset 
/ C (with N = {0, 1, 2, . . . }), and 9 i = d 1 ^ ■ ■ ■ d% . Let us denote by W u with I G 
L = {1, . . . , 2^^!}, the elements of the realization of the Weyl group Wn generated by the 
operators Ily and IL ■ The action of the total antisymmetrisation operator A over a factored 
state tp = <p(x) \s), with ip e C°°(R N ) and |s) € 5, is given by 

A^ = ^e i (W / /^)W / /|s), (21) 

where e; = ±1 is the parity of the total number of generators ELy, 11^ in any decomposition 
of Wi. By hypothesis, 

B(A^)= eifiW&iWw) -BiiW^s)) =0. (22) 

Applying the latter equation to a family of functions {^(x)} - eIxL satisfying the condition 



ijeixL 

r -\i£l,l£L 

det 



^(Wi^Ol^'^^O, (23) 
J jeixL 



we obtain 

B t (Wi\s)) = , for all i£l,l€L. 

In particular (taking / G L so that W/ is the identity) S,|s) =0 for all i € I. Since |s) £ <S 
is arbitrary, it follows that i?i = for all i £ I, and hence B vanishes identically. □ 

Lemma 2. If B e Q <S> & commutes with A then (AB)* = A* B* for all Ae® ® 



Proof. Using Eq. (19) and the hypothesis repeatedly we obtain: 

(AB)*A = ABA = AAB = A* AB = A*BA = A*B*A . 

The statement follows from the previous lemma. □ 

We shall often make use of the following immediate consequence of Lemma |[ 

Lemma 3. If A, B e D <8 £ commute with A then [A, B]* = [A*,B*]. 

We shall now construct a complete family of commuting integrals of motion for the hy- 
perbolic Bn Sutherland spin Hamiltonian ([!]). The construction is based on the observation 
that the operator H in Eq. Jl6| ) can be expressed as 

H = -Y J Jf- (24) 

i 

By ([?]), the operators 

i 

commute with one another. In view of the previous lemma, it suffices to prove that I p 
commutes with the total antisymmetriser A for all p € N to conclude that the star operators 
Ip (p € N) form a commuting family of integrals of motion of H* = —I*. We shall in fact 
prove the following stronger result: 
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Lemma 4. The operators I p (p £ N) commute with the permutation and sign changing 
operators Kij and Ki. 

Proof. First of all, the elementary permutation i^^i+i commutes with I p for i — 1, . . . , N— 1. 
Indeed, 2£i,j+i commutes with J, for j ^ i, i + 1 by (g), while for j = j, i + 1 we have 

2p-l 
r=0 

2)9-1 

r=0 

Since an arbitrary permutation can be expressed as the product of elementary permutations, 
this shows that Kij commutes with I p for all i j. Secondly, the sign reversing operator 
Km commutes with I p for all p, since 

K N Ji = J l K Nl if i < N, 

while 

KnJn = —JnKnJn — 2(6 + b')J]y = J^Km ■ 
This implies that Ki commutes with I p for an arbitrary i = 1, . . . , N, since 
= K iN [K N K iN ,I p ] — K iN [K iN K i7 I p ] = [Ki, I p ] . 

□ 

The operators /* (p £ N) thus form an infinite commuting family. Moreover, by examining 
the terms of highest order in the partial derivatives one can easily conclude that the set 
{Ip}p=i is algebraically independent. We have thus proved the main result of this section: 

Theorem 1. The operators {I p } p= i form a complete family of commuting integrals of mo- 
tion for the hyperbolic Bn Sutherland spin Hamiltonian H* — —I*. 

We also note that the constants of motion I* (p £ N) commute with the total permutation 
and sign changing operators Iljj and IT. This is a consequence of Lemma |] and the following 
general fact: 

Lemma 5. If A £ 2) <S> & commutes with K^ (resp. Ki) then A* commutes with Iljj 
(resp. Hi). 

Proof. We can write 

A = D-yK-y > (26) 

where K~ i is a monomial in Kki and Ki, D 1 £ D, and T is a finite set such that {K 7 | 7 £ T} 
is linearly independent. By hypothesis 

A = KijAKij = (D 7 ) KijK^Kij (27) 

where Kij{D 1 ) is the image of D 1 under the natural action of Kij in 2). Comparing ( |27| ) 
with d26| ) we conclude that for each 7 £ T there exists 7' £ T such that K i jK 1 K i j = K 7 >, 
and 

K l3 (D 7 ) = D Y . 
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On the other hand we have 

no- A*n« = K a( D i) s^Sh = J2 K v(PJ K i' = E D 7' «y • 

7GT 7GT 7GT 

Since (7')' = 7, we have V = T, and therefore the right-hand side of the previous formula 
equals A*. The equality Ilj^lTIj = A* is established in a similar way. □ 



III. SPECTRUM OF THE SUTHERLAND SPIN MODEL OF B N TYPE 

We shall now study the algebraic sector of the spectrum of the -Bjv-type Sutherland spin 
model (|l|) . The starting point in our discussion is the invariance under H of the space TZ m 
for all m — 0, 1, ... , which is an immediate consequence of Eq. ( pi|) and the definition of 
the operators J, . We shall construct a basis of the iJ-invariant space lZ m with respect to 
which the matrix of H\n m is upper triangular, thereby obtaining an exact formula for the 
spectrum of this operator. 

To derive the spectrum of H* from that of H, we shall make use of the identity 

H*[A(<p\s))]=A[(H<p)\s)], (28) 

wherew £ C°°{R N ) and \s) G S. The latter identity, which is an immediate consequence of 
Eq. ([jl]) and Lemma |J, implies that the spaces 

M m = A(Km <8 S) , m = 0,l,..., (29) 

are invariant under H* . From the basis of lZ m triangularizing H\n m we shall construct 
a basis of M m with respect to which H*\_\4 m is also represented by an upper triangular 
matrix. In this way we shall determine the spectrum of the Sutherland spin model of Bn 

type ©. 

Let us start by computing the spectrum of the operator H . Following closely the approach 
of Ref. j2| for spin models of An type, we shall define a suitable partial ordering in the set 
of (scaled) exponential monomials 

/„(x) = exp [2 njXj) , n = (rii, . . . ,njv) , -™<n j; <m, (30) 

i 

spanning the subspaces lZ m . We shall then show that the operator H is represented by a 
triangular matrix in any partially ordered basis of lZ m . 



The partial ordering in the basis (30) is defined as follows. Given a multiindex n = 
(ill, . . . , ripf) s we define the nonnegative and nonincreasing multiindex [n] by 

[n] = (\n il \,...,\n iN \) , where |n n | > ■ ■ • > \n lN \ . (31) 

If n, rtl € [Z N ] are nonnegative and nonincreasing multiindices, we shall say that n -< n 1 if 
ni — n'i = ■ ■ ■ = rii-i — n' i _ 1 — and rii < n^. For two arbitrary multiindices n, n' £ Z^, 
by definition n -< n' if and only if [n] -< [n'] . The partial ordering -< in induces a partial 
ordering in the exponential monomial basis ([30|), namely /„ ~< f n i if and only if n -< n'. The 
action of the Weyl group on the basis ( |30| ) preserves this partial ordering, i.e., if /„ -< /„' 
then Wf n -< Wf n < for all W G W N . 

If n = (ni, . . . , njy) G [Z^] and s G Z, we shall use the following notation: 

#(s) = card{i : n t = s} , 
£(s) = min{i : rii = s } ? 

with £(s) = +00 if rii 7^ s for all i = 1, . . . , N. For instance, if n = (5, 2, 2, 1, 1, 1, 0) then 
#(1) = 3 and £(1) = 4. The computation of the spectrum of H is based on the following 
result: 
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Proposition 1. If n € [Z N ] is a nonnegative and nonincreasing multiindex, the following 
identity holds: 



Jifn — ^n,ifn + ^ C n i f n ' , 



(32) 



n'ez" 

n n 

where c™ i £ R and 

^ f2n i + 6 + 6' + 2a(iV + i + l-#(n i )-2^(ni)), > 

™ ~~ \-b - b' + 2a(i -N), n, = 0. 

Proof. After some algebra one readily obtains the following expression: 



(33) 



Jifn — fn 



2n t + b + b' + 2a(N - 1) - 2a 



i - 1 



Oij - 1 



Aj - 1 



Ct- — 1 ■ - — 1 



-1— 2n< 



ay - 1 



A, - 1 



(34) 



26 : 



Z,; - 1 



-1 26 ,^- 2 -+l 



where 



z,- = e 



2i + 1 



ay 



Consider, for instance, the first term in ay in Eq. (p4|). Since j < z, ry > n^. If ry = Ui this 
term vanishes. If n, > rii we have 



Ct- — 1 



— r r 



(35) 



where the last sum only appears if rij — rii > 1 . In this case we have < max{7y —r,ni + r} < 
rij for all r = 1, . . . , ry — — 1, so the multiindices of the monomials in the summation 
symbol in Eq. ( |35| ) satisfy 

(ni, ...,nj-r,...,m+r,..., n N ) -< n . 

It may be likewise verified that the multiindices n' of the monomials arising from the re- 
maining terms in Eq. ( p4[ ) either coincide with n or satisfy n' -< n. The value of A n! i given 
in Eq. (33) can be computed by evaluating the constant part of the expression in square 
brackets in the right-hand side of Eq. ([34]). For instance, the first term in ay in Eq. ( |3^ ) 
contributes the quantity — 2a(£(ni) — l) to X n ,i- □ 

Note that Eq. ( |32| ) does not hold if n does not belong to [Z w ] , so that Proposition [l] in 
general does not determine the spectrum of the restriction of Ji to lZ m - On the other hand, 
for an arbitrary multiindex n G 7L N we shall only need the following weaker result: 



Corollary 1. If n € 1 then 



Jifn 



InA fn' i 



(36) 



[«']=< M 



for some real constants 7^ 
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Proof. We have 

■h!n = J*Wf[ n ] , 

where W is any element of Wn such that /„ = W f[ n \ - Eq. (|6|) then follows from the 
previous proposition, the commutation relations (|^)-(|l2|) and the invariance of the partial 
ordering -< under the action of the Weyl group. □ 

The algebraic spectrum of H can be computed in closed form using the previous results, 
which imply the following proposition: 

Proposition 2. For all n G Z N the following identity holds: 

Hf n = - £ Af n] + ]T ci f n , , with cj'el. (37) 

n -< n 

Proof. Let W be any element of Wn such that f n — W f[ n ]- Since H = — Ji commutes with 
W by Lemma EL from (pi) we obtain 



» n'eZ N ,i n'eZ N ,i 

n' [n] n' -< [n] 

Eq. ( |37| ) follows immediately from the latter equation, Corollary [l] and the invariance of the 
partial ordering -< under the action of the Weyl group. □ 

Let B m = {fn(j) I j = 1) • • • ! (2™ + 1) W } be any exponential monomial basis of the linear 
space 1Z m partially ordered according to -<, i.e., such that if n(j) -< n(k) then j < k. The 
previous proposition implies that the matrix of the restriction of H to lZ m with respect to 
B m is upper triangular. The eigenvalues of this matrix are its diagonal elements 

E n = -J2rfn],i\ -m< nj <m, j = l,...,N. (38) 

i 

It should be noted, however, that the algebraic eigenfunctions of H must satisfy appro- 
priate boundary conditions that we shall now discuss. In the first place, since the potential 
of the -BAr-type spin Sutherland Hamiltonian (|l|) diverges on the hyperplanes Xi ± Xj = 0, 
1 < i < j < N, as (xi ± Xj)~ 2 , we must require that the eigenfunctions of H vanish faster 
than (xi ± Xj) 1 / 2 near these hyperplanes. This yields the conditions (cf. Eq. (|l4|)) 

a,b>-. (39) 

Secondly, the eigenfunctions must be square-integrable on their domain, which (without loss 
of generality) shall be taken as the open set X C R w given by 

< x N < ■ ■ ■ < xi . (40) 

The algebraic eigenfunctions lying in lZ m will satisfy this condition if and only if 

^{b + b')+a(N -l) + m<0. (41) 

The latter inequality implies that the number of algebraic levels of H is finite, since m 
cannot exceed the integer mi defined by 

m 1 = max jm G N ~(b + b') + a{N - 1) + m < oj . (42) 
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Note, in particular, that there are no algebraic eigenfunctions unless the parameters in the 
potential verify the inequality 

^(b + b') + a(N -1)<0. (43) 

^From now on, we shall work on the maximal //-invariant subspace !Z m] . Remark 1. We 
could also have considered algebraic eigenfunctions of the Hamiltonian (|l|) antisymmetric 
under permutations but even under sign reversals. On these eigenfunctions, the action of 
the operators Sjj and Si coincides with that of the operators —Ky and -fQ, respectively. 
Therefore Eq. (fTq) in the definition of the star mapping should be replaced by 

(DK ai ---K ar )* = (-1Y' D S ar ■ ■ ■ S ai , 

where r' is the number of permutation operators in the monomial K ai . . . K ar . As a con- 
sequence, the Hamiltonian (|l|) is the image under the new star mapping of the operator 
H(-b,-b'), with H(b,b') given by Eq. ©. It follows from Eqs. © and @ that H* 
possesses algebraic eigenfunctions of even parity if and only if 

a>\, 6<~, ~(b + b')+a(N -1)<0. (44) 

In particular, these inequalities and Eqs. ( |39| ) and (Ej^) imply that H* cannot have both odd 
and even parity algebraic eigenfunctions for any values of the parameters. 
Let us turn now to the algebraic spectrum of the -Bjv-type Sutherland spin model (Q) . First 
of all, if Eqs. (^) and ( fl3| ) are satisfied the wavefunctions in Ai mi are normalizable and 
well behaved near the singular hyperplanes. Secondly, if ip £ lZ mi is an eigenfunction 
of H with eigenvalue E and |s) S S is an arbitrary spin state, it follows from Eq. ( |2S| ) 
that A((p\s)) G A4 mi is either zero or an eigenfunction of H* with the same energy E. 
The algebraic spectrum of H* is thus a subset of the algebraic spectrum of H . Instead 
of studying the conditions under which A(<p|s)) does not vanish, in the next proposition 
we shall directly construct from B mi a basis of M mi with respect to which the matrix of 
H*\M m is upper triangular. We shall use the following notation: 



N - M ,At1 

mo = 2MTT ' (45) 

where x denotes the smallest integer greater than or equal to x G K. 

Proposition 3. The H* -invariant space A4 mi is nonzero if and only if mo < mi. If this 
condition holds, a basis of M mi consists of states of the form 

A(f n \ Sl ,...,s N )), (46) 

where n G [Z and |si, . . . , sjv) satisfy: 

(2M + 1, if 0<m<mi 

i) *{m)<\- f - (47) 

[M, if i%i=0; 

ii) Si > Sj , if rii = rij and i < j ; (48) 
Hi) Si > , if n t = . (49) 

Proof. In the first place, since AW = e(W)A for any element W of the realization of WV 
generated by IUj, TEi , the space M mi is spanned by states of the form A(/„|s)), where 
n G [Z w ] and \s) G S. Moreover, from the definition of the total antisymmetriser A it 
follows that a state of the form (El) with n G \L N ] vanishes if and only if either s$ = Sj 
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when rii = rij > and i ^ j, or s, = ±Sj when = Uj = and i ^ j, or Sj = when 
rii = 0. In particular, the condition (|47|) is necessary to ensure that the state ( fig) does not 
vanish. Since this condition cannot hold if n-y < m , and n-y < my for all states in Ai mil it 
follows that A4 mi is trivial if m\ < tuq. 

On the other hand, if mo < mi all the states (ff(f)-(fl9|) are nonzero, and it is immediate to 
show that they are also linearly independent. Moreover, any nonzero state of the form ( |46| ) 
with n £ [Z N ] can be written as 

A[W(f n \ Sl , . . . , s N ))] = e(W) A(/„| Sl , . . . , s N )) , 

where n and \sy, . . . , sjv) satisfy (]47|)-([49|), and W E Wn is an element of the stabilizer of 

fn- □ 

Corollary 2. If mo < mi, i/ie dimension of the H* -invariant space A4 mi is given by 

a- Ska \ (mx{2M + l)+W\ 

dim(M roi ) = ( N . (50) 

Proof. Indeed, from Eqs. (^)-(^) it easily follows that 

j* ( k a \ V- (M\(2M + 1\ (2M + 1\ /mi(2M+l) + M\ 



AT +--+iV mi =JV 



□ 



The algebraic spectrum of the hyperbolic Bn Sutherland spin model ([!]) follows directly 
from Proposition |3| 

Theorem 2. If m < mi, the algebraic energies of H* are given by 

k = -E a «,*> ( 51 ) 

i 

where n E [Z w ] satisfies the condition ( [47[ ) awe? A n ,i is given by Eq. (|33|). 

Proof. Let ^>„ )S = A(/„|si, . . . , sat}) be an element of the basis ([46|)-(|49|) of M mi . Using 
Eqs. ( |37| ) and ( p8| ) we easily obtain 

ff> n , s = -X]A^„, s + £ cS'A(/ n ,| Sl) ... )SiV )). (52) 

The state A(/ n /|si, . . . , sjv)) is proportional to a basis element of the form ifjn",s', with 
jn" = Wf n ' for some element Vl^ of Wn and n" -< n. Therefore the matrix of H*\M mi m 
the basis (f46f)-(|49"|) is also upper triangular, with diagonal elements given by ([H]). □ 

The algebraic ground state of the Hamiltonian H* and its degeneracy d can be determined 
using Proposition ^ and the previous theorem: 

Proposition 4. The multiindex n yielding the algebraic ground state and the degeneracy of 
this state are given by 

/M x 

i) N < M: n = 0, d = 



N 



ii) N > M: r 2A/+1 2M+1 M 

n = (m , . . . ,m ,m - 1, • ■ ■ , m - 1, . . . , 1, . . . , 1, 0, . . . , ) , (53) 

where r = N- (m - 1)(2M + 1) - M. 
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Proof. Note, first of all, that from the definition (Eq) of rn it follows that 1 < r < 2M + 1. 
Let n S [Z^l be a nonnegative and nonincreasing multiindex satisfying nx < mi and 
Eq. ([47]) . The contribution to the algebraic energy (|5|) of all the terms A 2 i such that rij is 
equal to a fixed value k € {nj}jLi can be easily evaluated in closed form. Indeed, denoting 
for brevity 

»o =*(*:), *i =^A) + #(*)-!. 

1 fc 

a = (6 + b' + 2a(N - 1)) , a k = a , 

2a a 

from (|33|) we have, for k > 0, 



- X! A L = -4a 2 - to - »'i - "fc + I) 2 

i=i i=io 

= -4a 2 #(/c) [a 2 , + (i +«i - 2)a fc + | + i h + if) - ^(7i + 5ii) + 1 . 

(54) 

Similarly, the contribution to the algebraic energy of all the terms in Eq. (^) with k — 

-]TA 2 , = -4a 2 X> + «-l) 2 (55) 

i=i i=i 

is easily seen to equal the right-hand side of Eq. (|54|), since ao = a. The derivative of the 
right-hand side of Eq. (^4]) with respect to fc, with i and i\ fixed, is given by 

4a#(k)(2a k + i + ii - 2) . (56) 

This is clearly positive, since iq +ii — 2 > and > a — ^ > on account of (ff2|). Hence 
the energy decreases if k decreases, i$ and i\ being fixed. It follows that any multiindex n 
corresponding to the minimum value of the algebraic energy must be of the form 

n = (m, . . . , m, m — 1, . . . , m — 1, . . . , e, . . . , e) , (57) 

where e = or e = 1 and mo < m < mi . 

Let k be an integer in the range 1 to m. We shall consider next the change in the algebraic 
energy associated to the multiindex ( ]57| ) when #(fc) decreases by 1, while #(fc — 1) increases 
by 1 (including the case in which k = e = 1 and therefore #(fc — 1) = #(0) = 0). Suppose, 
for instance, that k > 2. Denoting i 2 — £(k — 1) + #(fc — 1) — 1, the change in the algebraic 
energy is given by 

4a 2 I - - i - i 1 ~ a k + 2) 2 - y)(i - h - i% - a k -i + l) 2 



»1 »2 



afe-i) 2 



(58) 



+ (i - io - h - oL k + l) 2 + {i-h-h- 

i—io i— 

= -4(1 + 2a(a fe + h - 1)) < -4(1 + 2aa k ) < . 

It may be similarly verified that when k = 1 and either e = or e = 1 the change in 
the algebraic energy is negative. This implies that the multiindex n yielding the algebraic 
ground state is of the form ( ^3|) if N > M, and zero otherwise. The degeneracy of the 
algebraic ground state then follows immediately from Proposition 0. □ 
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The algebraicground energy can be easily obtained from the previous proposition and 
Eqs. (|33] ) and (jSq). Indeed, denoting 

v = 2M + 1 , c = b + b' + 2m - a , 

the algebraic ground energy for even v (that is, half-integer M) is given by 

E^ e = - vm Q (4ml{av - 1) - 6cm - av - 2) + - (a 2 - 3c 2 )N + 2acN 2 ~ - a 2 N 3 , 
'3 3 3 

while for odd v (integer M), the algebraic ground energy reads 

Eq = Eq e + m (a + 2c + 2m (l — av)) . 

IV. THE Bjv-TYPE SUTHERLAND SPIN CHAIN 

In this section we shall introduce a quantum spin chain closely related to the hyperbolic 
Bn Sutherland spin model (0). We shall establish the integrability of this chain by explic- 
itly constructing a complete family of commuting integrals of motion associated with the 
integrals /* of the Hamiltonian (jl]). 

The starting point in this construction is the following expansion of the hyperbolic -B/v 
Sutherland spin Hamiltonian ([!]) in terms of the parameter a: 

H* = -J2 d ^ +aU * + a 2 U{x), (59) 

i 

where 

n* =J2[ sinlr 2 x ij 5 y + sinh ~ 2 x % «%] + Yl C 3 sinh_2 Xl ~ P' cosh ~ 2 x i) S i ( 60 ) 
fj(x) = ( sinh_2 x ij + sinh_2 <4) + (P 2 sinh" 2 x. ( - P' 2 cosh -2 ^) (61) 

and 

a a 

The Hamiltonian of the hyperbolic Sutherland spin chain of Bn type is by definition the 
operator H*°, where the superscript means that the coordinates X4 are replaced by the 
equilibrium points x® of the potential U, which satisfy the system 

BTI 

g>?,...,<) = 0, i = l,...,N. (62) 

A necessary condition for the system ([32]) to have a solution in the region Xj > 0, i = 
1, . . . ,N, is that P' 2 > (3 2 + 2(N — 1). In fact, there is strong numerical evidence that a 
solution exists if and only if \(3'\ > \j3\ + 2(N — 1). Note that this inequality corresponds to 
the condition ( [43] ) (when b > 0) or ( |44| ) (when 6 < 0) necessary for the existence of square- 
integrable algebraic eigenfunctions of the dynamical model, a fact certainly deserving further 
study. 

Let us define the operator J t E C°°(R N ) <g K by 

Ji = d Xi -aJi, i = l,...,N. (63) 
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We shall also denote 

2p = E^ 2P > P^ N - ( 64 ) 

Note that 2^ is the coefficient of a 2 in Ii — —H, and thus equals — U(x) by Eq. (|59|). We 
shall prove below that the operators {X*°}p =1 form a complete family of commuting integrals 
of motion for the Sutherland Bm spin chain Hamiltonian H*°. Let us begin by establishing 
the commutativity of the operators 2~*° for all p E N. In fact, the following stronger result 
holds: 

Proposition 5. The operators X* (p G N) form a commuting family. 

Proof. The proposition follows directly from the commutativity of the operators I* , taking 
into account that I* is the coefficient of a 2p in the expansion of I* in powers of a. □ 

We show next that X*° commutes with H*° for all p G N. Note that this is not a 
consequence of the previous proposition, since H*° is not proportional to Zj* = — [/(x°). 

Proposition 6. TTie operators X*° (p G N) commute with the Bn Sutherland spin chain 
Hamiltonian 7i . 

Proof. From (|^) it follows that 

[H,Ji] = 0, i = l,...,N. 

Using Eqs. (^) and (^9|) in the previous identity and equating to zero the coefficient of a 2 in 
the resulting expression we obtain 

dU 

[n ' Jl] = ~d^' i = i,...,N. 



It follows that 

N 2p-l 

r ) 1 dx 



The operator C\ = —\H 1 U(~x)\ vanishes identically, since TL does not contain partial deriva- 
tives and t/(x) is a symmetric even function of x. Note, however, that C p need not vanish 
for p > 1. Expanding in powers of a the identities [H, A] = [I p , A] =0 we obtain 

[H,A] = [2p,A] = 0, peN. 

By Lemma ^| we have 

[n*,x;]=c* p , pen. (65) 

^From the symmetry of the function U(x) with respect to permutations and sign changes it 
follows that dU/dxi commutes with Kjk and Kj for j. k ^ i, while 

OXi OXj OXi OXi 

^From these identities one can easily show that C* is of the form 

i 

The commutativity of H* with X*° follows from the latter equation, Eq. ( |65|) and the 
definition of the equilibrium points (|62|). □ 
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Finally, we prove the algebraic independence of the set {T*° }p =1 , thus establishing the 
integrability of the hyperbolic Sutherland spin chain of Bn type: 

Theorem 3. The operators {T*°}p =1 form a complete family of commuting integrals of 
motion for the Bn Sutherland spin chain Hamiltonian 7i*° . 

Proof. The set {Tp}p = i is algebraically independent, since the operators J® (i = 1, . . . , N) 
are linearly independent and commute with each other. The counterpart of Lemma || for the 
inverse of the star operator implies that the family {T^*}p =1 is also algebraically independent. 
The lemma now follows from the identity I°* =2*°. □ 

Note also that the constants of motion T*° (p € N) commute with the total permutation 
and sign reversing operators Iljj and IT . This follows from the identities [I*,ITj] = [/p,IT] = 
by taking the coefficient of a 2p . 
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